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PURPOSE:  The  purpose  of  this  Coastal  and  Hydraulics  Engineering  Technical  Note  (CHETN)  is  to 
present  results  of  initial  research  on  an  investigation  into  the  forcing  parameters  of  wave  runup  on  an 
intermediate,  barred  beach  using  innovative  observational  techniques.  A  brief  overview  of  the 
technology  and  processing  used  to  extract  wave  data  from  the  terrestrial  laser  scanner  is  presented,  and 
relationships  between  runup,  foreshore  change,  and  wave  characteristics  at  the  base  of  the  foreshore  are 
investigated.  In  addition,  discussion  is  provided  concerning  the  implications  and  causes  of  foreshore 
morphology  evolution  and  observed  tidal  oscillations  in  wave  observations  at  the  base  of  the  foreshore. 
A  new  equation  for  runup  is  also  presented.  It  is  based  on  a  mix  of  wave  parameters  measured  at  the 
base  of  the  foreshore  and  from  offshore.  This  new  equation  explains  73%  of  the  variance  in  wave  runup 
observations.  The  results  of  this  work  suggest  that  knowledge  of  both  the  surf-zone  and  beach 
morphology  prior  to  an  approaching  storm  is  critical  for  accurate  predictions  of  runup  inundation  and 
erosion. 

INTRODUCTION:  Accurately  predicting  beach  evolution  and  coastal  inundation  during  storms 
requires  models  that  correctly  predict  wave  runup  and  swash-zone  processes,  the  principle  drivers  of 
sediment  exchange  between  the  beach  and  surf-zone  (Masselink  and  Hughes  1998),  as  well  as  up-to- 
date  antecedent  data  to  initialize  conditions.  Recent  research  quantified  a  relationship  between  spatial 
variations  in  the  amount  of  wave-breaking  within  the  inner  surf-zone  and  the  volume  of  erosion  of  the 
sub-aerial  beach  during  a  nor’Easter  along  10  km  of  the  northern  Outer  Banks  (Brodie  and  McNinch 
2011).  Specifically,  regions  with  increased  wave  breaking  in  the  inner-surf  zone  (due  to  the  presence  of 
an  inner  bar)  were  found  to  correlate  with  regions  of  decreased  erosion,  perhaps  due  to  reduced  wave 
energy  impinging  on  the  beach.  In  this  technical  note,  results  are  presented  which  suggest  the  correlation 
may  be  explained  by  the  effect  that  wave  dissipation  over  variable  morphology  in  the  surf-zone  has  on 
the  resultant  energy  available  for  wave  runup,  and  therefore  inundation  and  erosion  of  the  beach.  In 
order  to  investigate  the  physical  drivers  of  wave  runup  in  the  field  during  storm  conditions,  an 
automated  terrestrial  laser  scanner  was  mounted  above  the  dune  (Figure  1A)  to  continuously  measure 
the  free  surface  profile  of  the  beach  and  water  surface.  This  unique  data  collection  technique  enabled 
simultaneous  collection  of  data  on  foreshore  evolution,  runup,  and  inner-surf-zone  wave  and  water  level 
characteristics  along  a  narrow  cross-shore  transect  (Figure  IB,  Brodie  et  al.  2010);  all  of  the  parameters 
forced  wave-driven  runup. 

Observations  of  the  beach  during  high-energy  conditions  and  post-storm  wrack-lines  have  long 
suggested  that  maximum  water  levels  during  storms  are  often  much  higher  than  measured  water  levels 
in  the  surf-zone  and  can  vary  significantly  alongshore.  A  growing  body  of  literature  suggests  that 
comparing  maximum  and  mean  water  levels  reached  during  storms  with  variations  in  beach  and  dune 
morphology  will  therefore  be  a  good  predictor  of  areas  of  the  coast  that  are  particularly  susceptible  to 
coastal  erosion  and  damage  during  extreme  storms  (Sallenger,  Jr.  2000;  Stockdon  et  al.  2007;  Plant  et  al. 
2010).  While  these  storm-impact  models  perform  better  than  random  models  (54%  accuracy  compared 
to  33%),  there  is  clearly  room  for  improvement  (Stockdon  et  al.  2007). 

Total  water  levels  reached  during  storms  are  calculated  as  the  sum  of  astronomical  tide,  surge,  wave- 
driven  setup,  and  wave  runup.  Wave-driven  setup  is  defined  as  the  super-elevation  of  the  mean  water 
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level  due  to  cross-shore  gradients  in  radiation  stress  induced  by  wave  breaking  (Longuet-Higgins  and 
Stewart  1964).  It  is  often  measured  at  the  shoreline  by  averaging  the  elevation  of  the  swash  over  a  given 
time  period,  (qr)  (Stockdon  et  al.  2006).  Wave  runup  is  the  time-varying  elevation  of  the  most 
shoreward  swash  excursion  on  the  foreshore,  with  maximum  runup  often  defined  as  the  2%  exceedence 
elevation  of  runup  swash  events  over  a  given  time  period  (Stockdon  et  al  2006;  Holman  1986). 

Relationships  between  runup  and  setup  and  important  forcing  variables,  such  as  beach  foreshore  slope, 
offshore  wave  height  and  period  (wavelength),  and  surf-zone  morphology,  have  been  the  focus  of  many 
studies,  and  the  relationships  found  often  varied  depending  on  the  type  of  beach  studied  (Holland  et  al. 
1995;  Raubenheimer  and  Guza  1996;  Guza  and  Thornton  1981;  Holman  and  Sallenger,  Jr.  1985; 
Nielsen  and  Hanslow  1991;  Ruessink  et  al.  1998;  Ruggiero  2004).  Stockdon  et  al.  (2006)  provide  a 
review  of  the  aforementioned  relationships  and  propose  an  empirical  relationship  between  the  2% 
exceedence  elevation  of  runup  (Ri)  and  foreshore  beach  slope  {[if),  local  wave  height  (e.g.,  10  m  water 
depth)  reverse-shoaled  to  its  deepwater  equivalent  (Ho),  and  deep-water  wavelength  (Lo): 


Eq.  1 

While  Equation  1  is  easy  to  calculate  for  given  boundary  conditions  ([if  Ho,  Lo),  accurately  defining  [if 
may  be  difficult,  particularly  on  intermediate  beaches  where  profiles  often  have  cross-shore  variations  in 
slope  and  can  change  on  the  timescale  of  hours  (Holland  and  Puleo  2001).  Though  Equation  1  is  used 
extensively  in  the  field,  a  surprising  amount  of  scatter  remains  unexplained  in  runup  data  sets.  Melby 
(2012)  also  reviewed  a  number  of  empirical  and  numerical  modeling  approaches  for  predicting  wave 
runup.  He  found  that  the  numerical  model  CSHORE  (Kobayashi  2009)  provided  similar  skill  to  the 
Stockdon  relationship  but  provided  significantly  poorer  results  on  dissipative  beaches.  For  both 
approaches,  he  noted  that  outdated,  forcing  morphological  data  may  have  skewed  results. 

Total  energy  in  the  swash  is  dependent  upon  contributions  from  both  the  incident  swell  band  (defined 
here  as  frequencies  (f),  0.05Hz  <  f  <0.25Hz)  and  infragravity  band  waves  (f  <  0.05Hz).  On  intermediate 
and  reflective  beaches,  the  incident  band  is  generally  more  dominant  in  swash  spectra  (Raubenheimer 
and  Guza  1996;  Holland  and  Holman  1993),  whereas  on  dissipative  beaches  or  during  storm  conditions, 
the  infragravity  band  tends  to  dominate  swash  variance  (Ruessink  et  al.  1998;  Ruggiero  2004).  Recent 
work  by  Guedes  et  al.  (2011)  suggests  that  the  tide  level  may  have  an  important  effect  on  swash 
statistics  (including  Rf),  particularly  at  intermediate  beaches  through  the  modulation  of  incident  band 
wave  breaking  over  sandbars.  This  dependence  on  conditions  within  the  inner  surf-zone  suggests  that 
wave  parameters  observed  outside  the  surf-zone,  such  as  those  used  in  Equation  1,  may  not  be 
universally  appropriate  for  runup  predictions.  Guedes  et  al.  (2011)  lacked  wave  data  within  the  surf- 
zone  and  consistent  beach  profile  observations  to  determine  if  the  tidal  oscillations  were  due  to  changing 
beach  slopes,  changing  wave  dissipation,  or  a  combination  of  both. 

The  objectives  of  this  study  were  to  further  test  relationships  between  runup,  morphology,  and  forcing 
conditions  in  the  field  using  a  terrestrial  laser  scanner  on  an  intermediate  beach  at  the  U.S.  Army  Corps 
of  Engineers  (USACE)  Field  Research  Facility  (FRF)  pier  in  Duck,  NC.  Data  from  the  terrestrial  laser 
scanner  are  unique  in  that  they  enable  simultaneous  observations  of  both  runup  and  important  forcing 
variables,  such  as  foreshore  slope  and  wave  conditions,  at  the  base  of  the  foreshore.  Similar  to  Guedes  et 
al.  (2011),  a  strong  tidal  dependence  in  de-tided  runup  statistics  was  observed,  which  directly  correlates 
with  tidal  variations  in  incident  band  energy  at  the  base  of  the  foreshore.  This  correlation  was  further 
altered  depending  on  beach  morphology  and  foreshore  slope,  which  was  observed  to  change  hourly  due 
to  profile  evolution  and  water  level  position. 
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METHODOLOGY:  A  description  of  the  data  collection  methodologies,  modeling  approach,  and 
experiment  conditions  are  presented  in  the  following  sections.  Some  brief  background  is  presented  on 
prior  approaches  for  measuring  wave  runup  to  highlight  the  advances  of  the  methodology  used  in  this 
study. 

Runup  Data  Collection.  Previous  studies  aimed  at  measuring  wave  runup  and  swash  zone  water 
levels  were  restricted  to  analyzing  water-elevation  time  series  of  (1)  the  shoreward-most  swash 
excursion  using  video  imaging  or  near-bed  resistance  wires,  (e.g.,  Holland  et  al.  1995;  Guza  and 
Thornton  1982),  or  (2)  the  free  water  surface  position  at  a  particular  location  on  the  foreshore  using 
pressure  sensors  (Raubenheimer  and  Guza  1996).  Recently,  approaches  using  lidar  technology  have 
been  developed,  which  allow  for  simultaneous,  high-frequency  observations  of  the  complete  free- 
surface  and  beach/swash  profile  position  (Brodie  et  al.  2010;  Blenkinsopp  et  al.  2010). 

The  basic  premise  of  lidar  technology  is  that  the  scanner  transmits  a  pulsed  light  beam  at  a  known  angle 
in  the  direction  of  a  target  and  records  the  time  of  travel  between  the  incident  and  reflected  pulses  to 
determine  the  distance  to  the  target,  and  when  combined  with  the  known  angle  of  the  pulse,  calculates 
elevation.  To  scan  along  a  narrow  vertical  profile,  termed  “line-scan,”  the  laser  beam  is  vertically 
deflected  by  a  polygon  of  spinning  mirrors  that  either  continuously  rotates  or  oscillates  up  and  down 
depending  upon  the  desired  scan  rate.  The  laser  can  also  provide  complete  360°  horizontal  coverage  by 
rotating  the  optical  head,  termed  “frame-scan.” 

For  this  experiment,  a  Riegl  terrestrial  laser  scanner  (LMS-z390i  1550pm  laser)  was  mounted  on  a  4  m 
guyed  scaffolding  tower  on  the  crest  of  the  primary  dune  roughly  400  m  north  of  the  FRF  pier  in  line 
with  the  FRF’s  cross-shore  array  of  wave  gauges  (Hanson  et  al.  2010),  Figure  1A.  The  instrument  is 
positioned  approximately  8  m  above  mean  sea  level.  The  laser  is  automated  to  complete  an  hourly  20 
minute  cross-shore  line-scan  at  2.5Hz  and  a  0.05°  angular  resolution,  followed  by  a  180°  frame  scan. 
Twenty  minute  collections  were  chosen  to  be  similar  to  the  commonly  used  Argus  eighteen  minute 
collection  length  (Lippmann  and  Holman  1989);  however,  this  collection  length  may  have  prevented 
inclusion  of  longer  period  infragravity  waves  in  this  analysis.  Since  the  laser’s  vertical  field  of  view  is 
+/-  40°  from  horizontal,  it  is  mounted  at  a  30°  angle  toward  the  ocean  to  enable  viewing  of  the  dune  toe. 
Though  angles  of  incidence  to  the  water  are  relatively  low  from  the  laser’s  position  on  the  dune,  the 
laser  reflects  well  off  of  the  water  surface  when  foam  is  present  (blue  dots,  Figure  IB).  Maximum  range 
of  measurement  depends  upon  the  amount  of  breaking  and  foam  present  in  the  surf-zone  at  any  given 
time,  but  rarely  exceeds  150  m  for  this  laser  scanner.  Drawbacks  to  the  laser  include  attenuation  in 
heavy  rain  or  fog,  and  shadowing  of  the  backside  of  incoming  waves  due  to  the  low  grazing  angle, 
which  therefore  requires  interpolation  to  generate  the  full  free-surface  profile.  In  addition,  returns  can  be 
sparser  on  the  foreshore  immediately  following  rundown  when  a  thin  film  of  water  is  left  on  the  beach 
surface.  Field  precision  and  accuracy  are  found  to  be  roughly  twice  manufacturer  specifications  at  0.015 
m. 

To  extract  foreshore  slope,  runup,  and  wave  statistics,  the  remainder  of  this  paper  focuses  mainly  on 
hourly  line  scans.  Data  are  first  transformed  from  the  scanner’s  own  coordinate  system  (angle,  range 
space)  to  rectified  Cartesian  coordinates  (local  FRF  coordinates  for  the  horizontal,  and  the  North 
American  Vertical  Datum  of  1988  (NAVD88)  coordinates  for  the  vertical),  using  rotation  matrices 
determined  from  scans  of  reflectors  with  precisely  known  coordinates.  Profile  data  are  then  interpolated 
to  10-cm  resolution  in  the  cross-shore  at  each  time-step.  To  separate  foreshore  position  from  the  free- 
surface  position  at  each  time-step,  a  1.2  sec  through  30  sec  running  minimum  foreshore  within  each  20 
minute  time  series  is  calculated,  and  data  points  are  classified  as  water  if  they  are  >0.015  m  above  the 
running  minimum  foreshore  (Figure  2).  A  running  minimum  foreshore  is  used  due  to  foreshore 
elevation  changes  that  occur  on  the  timescale  of  incident  and  infragravity  band  swash  (Howd  and 
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Holman  1987).  The  runup  time  series  is  therefore  defined  using  the  most  landward  water  point  at  each 
time-step  (Figure  2). 


Cross  Shore  Distance  (m) 


Figure  1 .  The  laser  is  shown  mounted  in  the  dune  in  panel  A,  and  an  example  of  laser  returns  off  the  water 
surface  (blue  dots)  over  twenty  minutes  is  shown  in  panel  B. 


Time  (s) 

Figure  2.  An  example  laser  time  stack  with  digitized  runup  time  series.  Colors  show  elevation  difference  from 
minimum  foreshore.  Points  with  elevation  differences  <  0.015  m  are  colored  black.  Laser  data  offshore  of  the 
minimum  foreshore  are  colored  red  for  less  confusion.  The  runup  line  is  shown  in  magenta. 

Surge  and  tide  are  removed  from  runup  time  series  by  subtracting  elevations  recorded  by  the  National 
Oceanic  and  Atmospheric  Administration  (NOAA)  tide  gauge  on  the  end  of  the  FRF  pier  where  wave 
setup  is  not  present.  The  resultant  wave-driven  swash  is  then  averaged  to  find  rjr.  Spectra  are  also 
calculated  on  the  wave-driven  runup  time  series  to  find  significant  infragravity  band  swash,  Ssig,  and 
significant  incident  band  swash,  Ssin.  In  addition,  observed  R2  is  calculated  for  each  runup  time  series 
based  on  the  method  described  by  Stockdon  et  al.  2006.  Wave  setup  ( rjw )  and  spectra  are  also  calculated 
at  the  base  of  the  foreshore,  defined  as  the  position  of  minimum  rundown  during  the  20  minute  line- 
scan,  for  each  collection.  These  wave  spectra  are  analyzed  to  calculate  significant  infragravity  band 
wave  height,  Hstg,  and  significant  incident  band  wave  height,  Hsm,  and  describe  the  hydrodynamic 
conditions  immediately  seaward  of  the  swash  zone.  As  described  in  Stockdon  et  al.  (2006),  /?/  is 
calculated  for  each  time-step  by  fitting  a  linear  trend  to  the  mean  beach  topography  during  each  20 
minute  line-scan  between  the  location  of  the  mean-swash  +/-  two  swash  standard  deviations  in  the  cross¬ 
shore.  Cross-shore  variations  in  beach  slope  were  calculated  by  taking  the  first  derivative  down  the  20 
minute  averaged  beach  profile  and  used  to  investigate  beach  concavity  at  each  time-step. 
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Empirical  Modeling  Approach.  For  comparison  to  observed  data.  Equation  1  was  calculated  for 
each  time-step  twice — once  using  the  observed  beach  slope  and  once  using  a  “mean  beach  slope”  for  the 
entire  experiment,  to  determine  the  sensitivity  of  the  equation  to  the  correct  beach  slope  input.  Wave 
parameters  ( Ho,Lo )  were  determined  by  reverse-shoaling  data  from  the  FRF’s  11m  Acoustic  Wave  and 
Current  (AW AC)  profiler  to  deep  water  values.  Local  water  levels  (tide  and  surge)  were  measured  at 
NOAA’s  tide  station  on  the  end  of  the  FRF  pier.  Degrees  of  freedom  for  correlation  analyses  were 
calculated  by  dividing  the  total  number  of  samples  by  12  to  conservatively  address  any  issues  of 
temporal  autocorrelation  within  the  time  series. 

Experiment  Conditions.  Analyzed  data  spanned  a  period  of  roughly  650  hours  of  semi-continuous 
laser  observations  from  26  Aug  through  29  Sep  2011,  including  the  landfall  of  Hurricane  Irene  and 
offshore  swell  from  Hurricane  Katia  and  Tropical  Storm  Maria,  and  from  26  Oct  through  13  Nov  2011, 
during  the  arrival  of  multiple  nor’easters.  Wave  conditions  during  the  experiment  as  recorded  by  the  1 1 
m  AWAC  at  the  FRF  are  shown  in  Figure  3,  with  red  shading  indicating  laser  collection  times.  Though 
time  gaps  exist  in  the  data  set,  for  ease  of  plotting,  most  graphs  herein  will  show  collection  number 
rather  than  time  on  the  x-axis,  and  vertical  magenta  bars  will  represent  breaks  in  time  (Figure  3B). 


Figure  3.  Wave  height  at  the  1 1  m  AWAC  (blue)  and  from  the  laser  (red)  is  shown  in  the  top  panel  with  time 
on  the  x-axis,  and  again  in  the  bottom  panel  with  laser  collection  number. 

RESULTS:  Experiment  results  are  presented  in  the  following  sections.  Morphology  changes  to  the 
beach  profile  are  discussed  in  the  first  section,  with  a  specific  concentration  on  variations  in  observed 
foreshore  slope.  Wave  observations  from  offshore  and  at  the  base  of  the  foreshore  are  discussed  next, 
followed  by  modeled  and  observed  swash  data. 

Foreshore  Slope.  Changes  in  foreshore  slope  during  the  course  of  the  experiment  are  analyzed  in 
Figure  4.  It  is  important  to  note  that  the  cross-shore  area  over  which  foreshore  slope  is  calculated  can 
change  depending  upon  the  location  of  the  swash  zone,  and  therefore  changes  in  foreshore  slope  can  be 
due  to  one  of  two  processes  (or  a  combination  of  both):  (1)  changes  in  tide  level  allowing  the  swash  to 
cover  a  steeper  or  flatter  portion  of  the  beach  or  (2)  bed-level  changes  in  the  foreshore  due  to  erosion 
and  accretion. 

Mean  foreshore  slope  during  the  experiment  was  0.1  m/m,  with  values  ranging  from  0.02  m/m  up  to 
0.22  m/m  (or  a  range  of  11.3°).  A  strong  positive  correlation  is  measured  between  water  level  and 
foreshore  slope  during  the  first  275  laser  collections  (blue  line,  Figure  4A),  in  which  steeper  slopes  are 
observed  at  higher  tides,  and  flatter  slopes  are  observed  at  lower  tides.  During  collections  275  through 
350,  little  relationship  between  foreshore  slope  and  water  level  is  observed.  Beginning  at  collection  350 
and  continuing  to  collection  450,  foreshore  slope  variation  is  out  of  phase  with  the  tide,  such  that  flatter 
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slopes  are  observed  at  high  tide  and  steeper  slopes  at  low  tide,  though  the  variations  in  slope  are  less. 
From  collection  450  through  550  and  600  through  the  end  of  the  experiment,  the  relationship  returns  to 
in-phase. 


A.  Foreshore  Slope  (blue)  and  Tide  (green) 
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Figure  4.  Analysis  of  changes  in  foreshore  slope  during  the  experiment.  Foreshore  slope  and  tide  elevation 
are  shown  in  panel  A.  In  panel  B,  a  time  stack  of  cross-shore  variations  in  slope  is  shown  with  warmer  colors 
representing  steeper  slopes  and  cooler  colors  representing  flatter  slopes.  In  panels  C-E,  example  profiles 

from  low  and  high  tide  are  shown. 


To  better  illustrate  these  variations  in  foreshore  slope,  hourly  cross-shore  beach  slopes  (colors)  are 
plotted  in  a  time  stack  in  Figure  4B,  and  example  high  and  low  tide  profiles  are  plotted  in  Figure  4C-E. 
During  the  periods  of  strong  in-phase  correlations  with  water  level,  the  beach  was  both  (1)  convex  in 
shape — steeper  slopes  (warm  colors)  higher  on  the  beach,  and  flatter  slopes  (cool  colors)  lower  on  the 
beach  (Figure  4B),  and  (2)  became  even  more  convex  at  high  tide  due  to  erosion  of  the  upper  foreshore 
with  the  rising  tide  (red  line,  Figure  AC).  When  a  lack  of  relationship  with  water  level  was  observed,  the 
foreshore  was  more  linear;  while  accretion  did  occur  with  the  rising  tide,  magnitudes  were  smaller  and 
more  evenly  distributed  across  the  profile  (Figure  4D).  Finally,  during  periods  of  weak  out-of-phase 
relationship  between  the  tide  and  foreshore  slope,  accretion  occurred  with  the  rising  tide  on  the  middle 
and  lower  foreshore.  This  caused  the  slightly  concave  low  tide  profde  (blue  line,  Figure  4E)  to  become 
slightly  convex  at  high  tide  (red  line,  Figure  4E). 

Observed  Wave  Data.  Offshore  wave  forcing  parameters  (sqrt(HoLo))  are  plotted  in  Figure  5A. 
Waves  during  Hurricane  Irene  (collections  1-50)  and  from  the  offshore  passage  of  Hurricane  Katia 
(collections  50-80)  as  well  as  during  the  nor’easters  at  the  end  of  the  experiment  (collections  600-650) 
were  characterized  by  high  values  of  sqrt(HoLo),  and  thus  strong  potential  for  high  runup  and  increased 
breaking  in  the  surf-zone. 
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Hourly  and  12  hour  averaged  measurements  of  de-tided  rjw  at  the  base  of  the  foreshore  (skinny  and  thick 
blue  lines,  respectively)  are  plotted  with  tide  (dotted  black  line)  in  Figure  5B.  Interestingly,  a  tidal  signal 
was  observed  in  de-tided  rjw  that  was  roughly  90°  out  of  phase  with  the  tide:  elevated  //„■  (setup)  was 
observed  at  mid-falling  tide,  whereas  lower  rjw  and,  in  some  cases,  negative  rjw  (setdown)  was  observed 
at  mid-rising  tide.  Averaging  over  these  tidal  variations  suggests  setup  at  the  base  of  the  foreshore  was 
greatest  during  the  nor’easters  at  the  end  of  the  experiment,  and  only  slightly  enhanced  during  Hurricane 
Irene  and  Katia  waves. 


A.  Offshore  Wave  Forcing:  sqrt(H0LQ) 
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Figure  5.  Offshore  wave  forcing  from  the  1 1  m  AWAC  (panel  A),  laser-observed  water  level  (panel  B),  and 
wave  height  characteristics  (panel  C)  at  the  base  of  the  foreshore  are  plotted  with  tides  during  the 

experiment. 

Significant  infragravity  band  wave  height  (blue  line)  and  significant  incident  band  wave  height  (red  line) 
at  the  base  of  the  foreshore  are  plotted  with  tides  (black  dotted  line)  in  Figure  5C.  Hsig  is  highest  during 
Hurricanes  Irene  and  Katia,  and  shows  slight  increases  during  Tropical  Storm  Maria  swell  and  the 
nor’easters  at  the  end  of  the  experiment,  similar  to  sqrt(HoLo),  as  expected.  A  strong  tidal  signal  is 
observed  in  Hsm  that  is  in  phase  with  the  tide  and  consistent  with  a  tidally-modulated,  saturated  wave 
field,  such  that  Hsm  is  larger  at  high  water  level  than  at  low  water  level  for  a  given  offshore  wave  height. 

Modeled  and  Observed  Swash  Statistics.  Runup  predictions  using  offshore  waves  and  foreshore 
slope  measurements  are  shown  in  Figure  6A,  along  with  phase-locked  modeled  (red  and  black  lines), 
observed  (blue  line),  R2,  and  water  level.  Wave-driven  runup  was  highest  during  Hurricanes  Irene  and 
Katia,  and  during  the  nor’easters  at  the  end  of  the  experiment.  A  strong  in-phase  tidal  signal  is  apparent 
in  de-tided,  observed,  R2  such  that  R2  is  higher  at  high  tide  than  low  tide  for  a  given  offshore  wave 
condition.  This  tidal  signal  is  only  apparent  in  the  model  data  when  the  observed  foreshore  slope  is 
included  (red  line),  and  therefore  during  periods  of  time  when  foreshore  slope  was  not  strongly 
correlated  with  water  level  (Figure  4A),  no  tidal  variation  was  predicted  by  Equation  1.  While  including 
observations  of  foreshore  slope  in  Equation  1  increased  the  correlation  between  predicted  and  observed 
data  from  an  r2  =  0.36  (p«0.001)  to  an  r2  =  0.56  (p«0.001),  however,  inclusion  of  slope  variation 
alone  does  not  result  in  sufficient  tidal  modulation  when  compared  with  measurements. 
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Figure  6.  Observed  and  modeled  swash  statistics  plotted  with  tides  during  the  experiment.  The  laser 

collection  number  is  shown  on  the  x-axis. 

The  results  for  r\R  are  similar  to  those  of  R2,  with  a  strong  in-phase  relationship  between  observed  data 
and  water  level  that  was  not  predicted  by  Equation  1 .  In  addition,  a  significant  correlation  between 
observed  data  and  tidal  elevation  exists  for  Sin,  with  the  exception  of  collection  numbers  250-450, 
where  the  relationship  mirrors  that  of  foreshore  slope  and  tides  (Figure  4A).  In  general,  there  seems  to 
be  little  tidal  dependence  in  Sig,  which  instead  correlates  well  with  sqrt(H0L0)  (r2= 0.65 ,  p«0.001). 

DISCUSSION:  The  above  results  are  discussed  in  the  following  sections  with  a  focus  on  understanding 
the  observed  tidal  variations  in  many  of  the  measured  variables.  The  results  suggest  that  increased  or 
decreased  breaking  over  the  sandbar  at  low  and  high  tides,  respectively,  may  modulate  energy  in  the 
incident  band,  affecting  runup  at  the  shoreline.  The  implications  of  this  sandbar  filtering  are  (1)  that  up- 
to-date  bathymetry  or  a  quantification  of  dissipation  across  the  surf-zone  in  addition  to  sub-aerial  beach 
morphology  is  necessary  to  produce  accurate  runup  predictions  that  capture  alongshore  variations  along 
the  coastline.  A  new  equation  for  predicting  wave  runup  is  proposed  that  includes  incident  band  wave- 
height  at  the  base  of  the  foreshore  with  beach  foreshore  slope  and  offshore  wave  parameters,  and  it 
explains  73%  of  the  variance  in  the  observations.  The  feasibility  of  using  this  equation  in  practice  is  also 
addressed. 

Foreshore  Slope.  Changes  in  both  erosion/accretion  and  the  cross-shore  position  of  the  swash  with 
the  rising  and  falling  tide,  as  well  as  changing  wave  conditions,  can  contribute  to  the  hourly  evolution  of 
foreshore  slope  on  intermediate  beaches.  In  this  experiment,  the  relationship  between  erosion  and 
accretion  patterns  and  the  rising  and  falling  tide  changed  depending  upon  the  shape  of  the  beach.  A 
concave  beach  that  featured  steep  slopes  on  the  upper  foreshore  and  flatter  slopes  on  the  lower  foreshore 
was  further  eroded  during  the  rising  tide,  creating  even  steeper  slopes  in  the  swash  zone.  In  contrast,  a 
slightly  convex  beach  that  featured  flatter,  upper  foreshore  slopes  and  steeper,  lower  foreshore  slopes 
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tended  to  accrete  with  the  rising  tide,  creating  even  flatter  slopes  within  the  high  tide  swash  zone.  Linear 
beaches  tended  to  accrete  and  erode  more  uniformly  with  the  rising  and  falling  tide. 

Though  profile  elevation  changes  were  only  on  the  order  of  +/-  0.15  m,  and  the  mean  cross-shore 
position  of  the  swash  only  varied  by  roughly  10m  from  high  tide  to  low  tide,  changes  in  slope  during  in- 
phase  collection  times  (convex  beach)  often  reached  4°,  a  change  of  -70°  of  the  mean  beach  slope. 
Surprisingly,  a  change  in  4°  of  slope  will  change  the  runup  predictions  made  by  Equation  1  by  47%, 
which  highlights  the  potential  errors  in  using  Equation  1  when  the  current  beach  foreshore  slope  is  not 
known. 


Wave  Climate  at  the  Base  of  the  Foreshore.  Significant  wave  height  in  the  incident  band  at  the 
base  of  the  foreshore  oscillated  strongly  with  the  tides,  such  that  higher  Hsm  were  observed  at  high-tides 
and  lower  Hsm  were  observed  at  low-tide  for  the  same  given  offshore  wave  conditions.  Video  imagery 
collected  by  the  Argus  cameras  at  the  FRF  was  used  to  investigate  wave  breaking  in  the  surf-zone  as  a 
possible  explanation  of  this.  An  example  set  of  imagery  is  shown  from  28  Sep  2011  at  high  tide  (Figure 
7A)  and  low  tide  (Figure  7B),  with  the  location  of  the  lidar  profile  in  magenta.  The  image  is  a  “bright” 
image,  which  saves  the  brightest  intensity  at  each  pixel  during  the  hourly  Argus  collection,  and  thus 
maps  regions  of  wave  breaking  (white  pixels)  in  the  surf-zone.  As  seen  in  panel  A,  at  high  tide  there  is  a 
small  amount  of  breaking  over  the  nearshore  sandbar,  and  a  wide,  energetic  swash  zone  is  observed  on 
the  beach.  In  contrast,  at  low  tide,  there  is  significantly  increased  breaking  over  the  sandbar  (due  to 
lower  water  levels),  and  a  more  narrow  swash  zone  can  be  seen  on  the  beach.  This  supports  wave  data 
from  the  lidar,  which  indicated  lower  wave  heights  at  the  base  of  the  foreshore  at  low  tide.  These 
observations  suggest  that  increased  breaking  of  incident  band  waves  at  low  tide  over  the  offshore 
sandbar  decreases  the  amount  of  incident  band  energy  reaching  the  shoreline,  supporting  one  hypothesis 
proposed  by  Guedes  et  al.  (2011).  Interestingly,  neither  Hsig  nor  SIG  correlates  with  water  level,  even 
though  increased  breaking  at  low  tide  should  promote  the  release  of  infragravity  band  waves  (Ruessink 
1998). 


Figure  7.  Argus’  brightest  pixel  images  for  a  high  tide  (A)  and  subsequent  low  tide  (B)  showing  variations  in 

wave  breaking  in  the  surf-zone. 


Though  ipv  at  the  base  of  the  foreshore  did  show  a  strong  tidal  oscillation,  it  was  roughly  90°  out  of 
phase  with  the  tide,  having  peaks  at  mid-falling  tide  and  valleys  at  mid-rising  tide.  This  suggests  that  the 
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rapid  increase  of  setup  at  the  shoreline  may  begin  slightly  farther  offshore  during  low  and  rising  tides, 
whereas  at  high  and  falling  tides,  wave  shoaling  in  the  shore  break  at  the  base  of  the  foreshore  may 
actually  create  a  slight  setdown  in  water  levels  before  the  rapid  rise  in  the  swash.  Changes  in  wave 
shape  and  the  development  and  movement  of  a  shoreline  “step”  feature  during  the  rising  and  falling 
tides  (Bauer  and  Allen  1995)  may  also  contribute  to  these  observed  patterns. 

Runup  Predictions.  Though  R2  predictions  from  Equation  1  were  improved  when  observed  changes 
in  beach  slope  were  included  on  an  hourly  basis,  it  still  only  explained  54%  of  the  variance  in  the 
observed  R2  data,  and  could  not  reproduce  the  observed  tidally  dominated  variation  in  runup  when  beach 
slope  was  not  in  phase  with  the  tide.  Examination  of  Equation  1  reveals  that  the  only  functional 
dependence  of  runup  on  water  level  is  through  the  inclusion  of  the  foreshore  slope  parameter,  because 
the  remaining  forcing  parameters  are  based  on  deep  water  wave  height.  However,  as  described  above, 
the  data  instead  suggest  that  the  incident  band  wave  conditions  at  the  base  of  the  foreshore  dominate 
runup  variations  on  an  intermediate  beach.  Specifically,  the  data  show  that  wave  conditions  at  the  base 
of  the  foreshore  can  change  with  water  level  independently  of  the  offshore  wave  conditions,  before 
interacting  with  the  foreshore,  due  to  shallow-depth  driven  wave  breaking. 


To  quantify  these  observations  and  their  effect  on  runup,  observed  R2  was  deconstructed  into  observed 
rjR,  Sin,  and  Sig  (similar  to  the  approach  of  Equation  1)  and  these  variables  compared  with  forcing 
parameters  measured  at  the  foreshore  (Bf  and  Hsm)  and  offshore  ( sqrt(HoLo )).  Specifically,  ///?  showed  a 
positive  correlation  with  Hsm  (r2  =  0.5,  p«0.001),  Sin,  showed  a  strong  positive  correlation  to  HsinBf  (r2 
=  0.73,  p«0.001),  and  Sig  showed  a  strong  positive  correlation  to  sqrt(HoLo)  (r2  =  0.66,  p«0.001). 
Combining  these  relationships  in  a  similar  way  to  Stockdon  et  al.  (2006),  yields  a  new  equation  for 
runup: 


R?  =  1.1  {  0.8 Hsin  +  ■ 


(25HsINBf)2+(0.06jHj^) 


Eq.  2 


A  scatter  plot  of  predicted  R2  with  Equation  2  and  observed  R2  is  shown  in  Figure  8.  Equation  2  explains 
73%  of  the  variance  in  observed  R2,  compared  to  Equation  1,  which  only  explained  56%  of  the  variance. 
In  addition,  removing  the  dependence  on  the  highly  variable  parameter,  Bf,  and  altering  the  coefficient 
before  Hstn  in  the  Sin  component  from  25  to  2.5  in  Equation  2,  explains  69%  of  the  variance  in  observed 
R2.  While  knowing  beach  slope  is  still  important  in  runup  predictions,  this  suggests  that  runup  can  be 
almost  as  well  predicted  by  only  knowing  the  morphology  of  the  surf-zone  and  its  effect  on  wave  energy 
reaching  the  beach.  In  fact,  Equation  2,  without  including  beach  slope,  explained  more  of  the  variance  in 
wave  runup  during  the  experiment  when  compared  with  the  commonly  used  Equation  1. 

Wave  height  in  the  incident  band  at  the  base  of  the  foreshore  is  rarely  known,  and  thus  Equation  2  may 
not  be  practical  for  predictive  use;  however,  it  better  represents  the  physical  drivers  of  runup, 
particularly  on  barred  intermediate  beaches.  A  more  readily  available  remotely  sensed  variable,  such  as 
a  measure  of  wave  breaking  and  bar  morphology  from  video  or  radar  imagery,  could  be  investigated  for 
use  as  a  scaling  factor  to  convert  offshore  wave  parameters  to  Hsm  at  the  base  of  the  foreshore,  making 
Equation  2  more  readily  applicable. 

Implications  for  Alongshore  Variations  in  Runup  and  Inundation.  This  study  did  not 
specifically  address  alongshore  variations  in  runup,  however,  the  temporal  variations  observed  can  be 
translated  to  provide  a  possible  explanation  of  the  correlation  between  sandbar  morphology  and  beach 
erosion  observed  by  Brodie  and  McNinch  (2011).  Specifically,  this  study  illustrated  how  temporal 
variations  in  wave  breaking  over  a  sandbar  during  the  course  of  a  tidal  cycle  in  combination  with  the 
shape  of  the  beach  could  create  temporal  variations  in  wave  runup  for  a  constant  offshore  wave  height. 
Similarly,  alongshore  variations  in  sandbar  configuration  will  dissipate  different  amounts  of  wave 


10 


ERDC/CHL  CHETN-l-92 
October  2016 


energy,  creating  alongshore  variations  in  the  amount  of  wave  energy  that  reaches  the  beach,  even  though 
the  offshore  wave  conditions  may  be  the  same  along  the  coastline.  Small-scale  alongshore  variations  in 
the  slope  and  curvature  of  the  beach  can  then  further  enhance  or  inhibit  runup,  creating  more  alongshore 
variation  in  runup.  A  schematic  summarizing  this  concept  is  presented  in  Figure  9,  where  panel  A 
depicts  a  multiple-barred  stretch  of  beach  with  increased  wave  dissipation  through  the  surf-zone  and  less 
runup.  Panel  B  depicts  steeper  surf-zone  morphology  with  less  breaking  over  a  single  bar  and  more 
energy  reaching  the  shoreline  and  interacting  with  a  steep  beach,  ultimately  leading  to  higher  runup  and 
increased  potential  for  erosion.  Therefore,  current  observations  or  model  predictions  of  the  antecedent 
morphology  are  paramount  to  being  able  to  correctly  predict  runup  inundation  and  refine  vulnerability 
assessments  immediately  prior  to  a  specific  storm  threat. 


r2  =  0.73 


Figure  8.  Example  GIS  window  showing  spatial  SSURGO  soils  data  on  top  of  FAO  soils  data,  along  with  the 

FAO  soil  attributes. 
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Figure  9.  Diagram  showing  effect  of  surf-zone  and  beach  morphology  on  wave  dissipation  and  runup. 

CONCLUSIONS:  A  dune-mounted  terrestrial  laser  scanner  was  used  to  make  hourly  observations  of 
beach  morphology,  swash,  and  inner  surf-zone  hydrodynamics  to  investigate  the  physical  drivers  of 
runup  on  an  intermediate  beach  under  a  variety  of  wave  conditions.  Data  revealed  a  strong  tidal  signal  in 
both,  the  wave  height  at  the  base  of  the  foreshore  and  the  resultant  de-tided  2%  exceedence  observations 
of  runup.  Specifically,  mean  and  incident  band  swash,  as  well  as  incident  band  wave  height  at  the  base 
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of  the  foreshore  were  observed  to  oscillate  in-phase  with  the  tide,  such  that  for  a  given  offshore  wave 
condition,  all  were  elevated  at  high  tide  relative  to  low  tide.  Coincident  video  imagery  suggested  that 
this  was  primarily  due  to  increased  breaking  of  incident  band  waves  at  low  tide  over  a  nearshore  bar, 
which  dissipated  wave  energy.  Based  on  these  relationships,  a  new  equation  was  developed  (Equation 
2),  and  shown  to  explain  20%  more  of  the  variance  in  observed  R2  when  compared  with  the  commonly 
used  Stockdon  et  al.  (2006)  equation  (Equation  1).  Equation  1  appeared  overly  sensitive  to  changes  in 
beach  slope,  which  was  observed  to  change  by  as  much  as  4°  from  high  tide  to  low  tide  and  11°  over  the 
course  of  the  entire  experiment.  Changes  in  foreshore  slope  were  due  to  small-scale  profde  adjustments 
with  tide  (erosion  and  accretion  on  the  order  of  +/-  0.15  m)  and  were  either  in  phase  with  the  tide 
(concave  profdes),  180°  out  of  phase  with  the  tide  (convex  profdes),  or  independent  of  the  tide  (linear 
profiles).  Interestingly,  even  though  beach  slope  was  highly  variable  throughout  the  experiment, 
removing  it  from  Equation  2  only  reduced  r2  values  by  6%,  suggesting  that  changes  to  the  wave  climate 
across  the  surf-zone  may  have  a  larger  effect  on  runup,  than  variations  in  beach  slope. 

These  data,  in  combination  with  prior  work  (Brodie  et  al.  2011),  suggest  that  knowledge  of  the 
morphology  of  the  profile  and  surf-zone  immediately  prior  to  an  approaching  event  is  critical  for 
producing  correct  predictions  of  wave-driven  runup  during  storms.  Furthermore,  extreme  or  long- 
duration  events  frequently  cause  significant  morphological  evolution  to  the  system,  which,  in  turn,  may 
alter  runup  predictions  during  the  storm.  Therefore,  effort  should  also  be  focused  on  using  numerical 
models,  such  as  CSHORE  (Kobayashi  2009),  to  model  the  potential  morphological  changes  to  the 
system,  so  that  runup  predictions  remain  accurate  throughout  the  course  of  an  event. 

ADDITIONAL  INFORMATION:  For  additional  information,  contact  Kate  L.  Brodie,  Coastal 
Observation  and  Analysis  Branch,  Coastal  Hydraulics  Laboratory  by  mail  at  1261  Duck  Rd,  Duck,  NC 
27949;  by  phone  at  252-261-6840  extension  233;  or  via  email:  Katherine.L.Brodie@usace.army.mil. 

This  CHETN  should  be  cited  as  follows: 

Brodie,  K.  L.  2016.  Importance  of  antecedent  beach  and  surf-zone  morphology  to  wave  runup 
predictions.  Coastal  and  Hydraulics  Engineering  Technical  Note.  ERDC/CHL  CHETN-I-92. 
Vicksburg,  MS:  U.S.  Army  Engineer  Research  and  Development  Center. 
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